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Abstract 

Using a sample of 2.59 x 10^ ^(25') decays collected by the CLEO-c detector, we present results 
of a search for the decay chain ip{2S) —>■ n^hc,hc —>■ n{n~^'K~)Tr^ ,n = 1,2,3. We observe no 
significant signals for n = 1 and n = 3 and set upper limits for the corresponding decay rates. 
First evidence for the decay he ir'^'!T~TT'^iT~iT^ is presented, and a product branching fraction of 
B{iIj{2S) -> he) X B{hc -> 2(7r+7r-)7r°) = 1.88l|];|^lg:^J x 10"^ is measured. This result implies 
that he hadrons and he Jije have comparable rates, in agreement with expectations. 

PACS numbers: 13.25.Gv,14.40.Gx 
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Although the field of charmonium spectroscopy is now thirty-five years old, data on the 
cc singlet state, the hc(^Pj), remains sparse. Two experiments have identified the he and 
accurately measured its mass. The CLEO [l|, [1] measurements were made using the decay 
chain ip{2S) vr^/ic, vr*^ 77, he — > 'jrjc, and identifying the he either by fully reconstructing 
the event using many different hadronic decay channels of the rjc, or by reconstructing the 
7r° and 7 in the decay chain and inferring the existence of the t]c. The E835 experiment ^ 
made scans of antiproton energy and observed the reaction pp ^ h^ 'yrjcfjc —>■ 77- 
The experiment also searched for the evidence of the he in the previously reported isospin- 
suppressed decay he n^J/ip but none was found. The he is also expected to decay 
directly to multi-hadron final states; however, such decays have yet to be observed. The he 
width into such states is expected to be, by coincidence, comparable to that of the radiative 
decays; Godfrey and Rosner [^] predict branching fractions of 38% for 777c decays and 57% 
ggg decays, with the remainder being 7(7(7. The he-, unlike the Xcj mesons, has negative 
G-parity, and thus its multi-pion decays are likely to involve an odd number of pions. Here 
we report the results of a search for the decays of the he into ?7,(7r"'"7r~)7r°, with n = 1, 2, 3. 

The data presented here were taken by the CLEO-c detector operating at the Cornell 
Electron Storage Ring (CESR) with e^e~ collisions at a center of mass energy corresponding 
to the ipi^S) mass of 3.686 GeV. The data correspond to an integrated luminosity of 56.3 
pb~^ and the total number of ip{2S) events, determined according to the method described 
in 0], is calculated as (2.59 ± 0.05) x 10^. Like the previous CLEO analyses, we search for 
the he mesons produced by the isospin-violating decay i'{2S) ii^he- 

Charged particles are detected in a cylindrical wire chamber system immersed in a 1.0 T 
axial magnetic field induced by a superconducting solenoid. The solid angle for detecting 
charged particles is 93% of 47r, and the resolution 0.6% at 1 GeV. To identify the pions, we 
measure the specific ionization, dE/dx, in the drift chamber and require that it be within 4 
standard deviations of that expected for a pion. Photons are detected using the Csl crystal 
calorimeter also inside the magnet coil, which has an energy resolution of 2.2% at 1 GeV and 
5% at 100 MeV. Photon candidates are required to have a lateral shower shape consistent 
with that expected for a photon and not to align with the projection of any charged particle 
into the calorimeter. We combine photon pairs to make vr" candidates, and kinematically 
constrain them to the known 7r° mass; combinations with a of l^ss than 10 for the one 
degree of freedom are retained for further analysis. 

For each decay mode, we combine the requisite number of charged pion candidates with 
one vr'' candidate to form an he candidate. These particles are kinematically constrained 
with the beamspot to form a primary event vertex. We then add a second 7r° candidate in 
the event, ensuring that no photon is used in both candidates, to make a ipi^S) candidate. 
This ip{2S) candidate is then kinematically constrained to the four-momentum of the beam, 
the energy of which is calculated using the known ip{2S) mass. The momentum is non-zero 
only due to the crossing angle (~ 3 mrad per beam) in CESR. To make our final selection, 
we require the ip{2S) candidate to have a of l^ss than 25 for the four degrees of freedom 
for this fit; this requirement rejects most background combinations. 

The kinematic fit produces an he mass resolution which is much improved over a di- 
rect measurement of M(n(7r+7r~)7r°) and slightly improved compared to a measurement of 
the missing mass using the measured parameters of the transition alone. To study the 
efficiency and resolutions, we generated Monte Carlo samples for each he decay using a 
GEANT-based detector simulation [s^. The decay products of the he were generated ac- 
cording to phase space. For each of the three muliti-pion decays sought, the MC studies 
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show that the he mass distribution is well-represented by a double Gaussian signal shape 
over a slowly varying background. For the n = 2 case, for example, the shape parameters 
are anarrow = 1-19 MeV, a^iide = 3.18 MeV, and Nnarrow / Ntotai = 0.643. The efficiencies 
are shown in Table I. 

TABLE I: For each he decay mode, the efficiency, the raw event yield with statistical uncertainties 
obtained from the fit to the data, and the product branching fraction Bi x B2, where Bi = 
B{ijj{2S) TT^hc), and B2 = B{hc — > n(7r+7r~)7r'^), including systematic uncertainties. Upper 
hmits are quoted at 90% confidence level, and include the effects of systematic errors as described 
in the text. 



Mode 


efficiency (%) 


Yield 


Bix B2X 10^ 


7r+7r~7r'^ 
2(7r+7r-)7r0 
3(7r+7r-)7r0 


27.0 
18.8 
11.5 


92^1 
35 ±26 


< 0.19 

n OQ+0.48+0.47N, 
(^i.OO_g 45_g 30; 

(1.2 ±0.9 ±0.3) (< 2.5) 



The final invariant mass distributions are shown in Figs. 1(a), 2(a) and 1(c). In the case 
of he TT'^Ti^Ti^ the events are dominated by ipi^S) J/ip, with the subsequent decay 

of the J/ijj into two charged particles. The J/ijj has a very large branching into yU+yU" and 
these events will, in general, pass all selection criteria and enter the plot (Fig. 1(a)). The 
most efficacious way of eliminating these events is to reject those events with 3.0 < M,r+7r- < 
3.2 GeV/c^. Figure 1(b) shows the plot after this cut has been made. Neither Figs. 1(a) or 
1(b) show any excess in the he region. These histograms are fit to a background function 
(second order polynomial for Fig. 1(a) and an ARGUS style background function [9| for 
Fig. 1(b)), and signal function of fixed mass and width; the he mass is taken from to find 
90% confidence level upper limits of < 94 and < 14 events respectively. 

Fig. 2(a) shows the invariant mass distribution for he 2{n^Tt~)T{^. It shows a distinct 
excess in the region of the he- The distribution is fit to an ARGUS style background function, 
plus a fioating mass signal with a fixed shape from the Monte Carlo studies. The measured 
peak mass is 3525.6 ± 0.5 MeV, which may be compared with the Particle Data Group 0] 
number of 3525.93 ±0.27 MeV and the more recent CLEO |i] measurement of 3525.28 ±0.22 
MeV. The yield is 92^22 events, and has a significance of AAa. We also analyzed a large 
sample of Monte Carlo events generated using the known decays of the ip{2S) and designed 
to mimic the real data sample. Those events where an he meson was generated are explicitly 
excluded. Figure 2(b) shows the 2(7r+7r^)7r° mass plot from the remaining events and, as 
expected, it shows no sign of an excess in the he region. This mass distribution falls slightly 
faster than the equivalent one in data, demonstrating the lack of complete knowledge of 
ipi^S) decays, but it can be well fit by an ARGUS type background function. 

Fig. 1(c) shows the mass distribution for he — * 3(7r+7r^)7r*^. It shows a small, but not 
statistically signficant, excess in the signal region. The fit shown uses the same fixed mass 
of the he and the measured yield is 35 ± 26 events, corresponding to a 90% confidence level 
upper limit of 70. 

We consider systematic uncertainties from many different sources, and these are listed 
for the 2(7r+7r^)7r° mode in Table II. We assign uncertainties of 0.3% and 2%, respectively, 
on the detection efficiency for each track and for each photon. The largest systematic 
uncertainty in the 2(7r"'"7r~)7r° mode is due to uncertainties in the fitting procedure. The 
fit is performed in small mass bins to minimize fiuctuations due to choice of binning, and 
has a per degree freedom of 242/235. Using a background function of a second order 
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Chebyshev polyomial gives higher yields but a less satisfactory fit. Fits are also performed 
over wider and narrower mass ranges and using higher order polynomial background shapes. 
The systematic uncertainty is calculated from observing the range of yields from different, 
reasonable, fitting procudures. The h^. is known to be relatively narrow and our Monte 
Carlo simulation assumed an intrinsic width of 0.9 MeV. We assign a systematic uncertainty 
based upon the variation of yield if this number was in the range 0-1.5 MeV. To evaluate 
the systematic uncertainty due to our knowledge of the resolution, we allowed for variations 
of up to 10% in the width of the resolution function. To account for possible substructure in 
the Stt decay products a series of Monte Carlo samples were generated where the vr mesons 
are the product of intermediate p mesons, and we look at the spread of different efficiencies 
calculated. 

To convert the yields to product branching fractions, we divide by the product of the 
number of ip{2S) events in the data sample and the efficiency from Table I. For evaluating the 
limits in the cases where there is no significant signal, we take the probability density function 
and convolve this with Gaussian systematic uncertainties. We then find the branching 
fraction that includes 90% of the total area. 

TABLE II: Systematic uncertainties for the 2(7r''"7r~)7r'^ mode. 

Source Uncertainty (%) 

Efficiency of tracks and photons 10% 
Background function and fitting range 11^% 
cut efficiency 4% 
Signal natural width 5% 
Signal resolution 8% 
Possible substructure 6% 
Possible decays to J/ip 't^yo 

jV(-0(2S)) 2% 

Tbt^^ tWo 



The product branching fraction, B{ip{2S) he) x B{hc 2(7r+7r )7r'^) is calculated to 
be (l.SStoillloil) X 10-5. We note that this is ^ 5% of B{iP{2S) -> he) x B{K ^ IVc) 0- 
Given the large number of different hadronic final states that are available for he decays, we 
can conclude that these hadronic states have a width the same order of magnitude as the 
radiative decays into the rjc- 
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FIG. 1: Invariant mass plots for (a) 7r"'"7r~7r° (b) tt+tt^tt^ with a J/V' veto (c) 3(7r"'"7r~)7r°. Fig 
1(a) is fit using a second order Chcbychcv polynomial shape background. Figs. 1(b), and 1(c) are 
fit using an ARGUS type background function. 
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FIG. 2: Invariant mass plots for 2(7r+7r )7r*' for (a) data, and (b) non-hc Monte Carlo events. In 
each case the background function is an ARGUS type function. 
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